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ABSTRACT 

This paper aims at explaining the two phases in the observed specific star formation 
rate (sSFR), namely the high (> 3/Gyr) values at z >2 and the smooth decrease 
since z = 2. In order to do this, we compare to observations the specific star formation 
rate evolution predicted by well calibrated models of chemical evolution for elliptical 
and spiral galaxies, using the additional constraints on the mean stellar ages of these 
galaxies (at a given mass). We can conclude that the two phases of the sSFR evo- 
lution across cosmic time are due to different populations of galaxies. At z > 2 the 
contribution comes from spheroids: the progenitors of present-day massive ellipticals 
(which feature the highest sSFR) as well as halos and bulges in spirals (which con- 
tribute with average and lower-than- average sSFR). In each single galaxy the sSFR 
decreases rapidly and the star formation stops in <1 Gyr. However the combination 
of different generations of ellipticals in formation might result in an apparent lack of 
strong evolution of the sSFR (averaged over a population) at high redshift. The z < 2 
decrease is due to the slow evolution of the gas fraction in discs, modulated by the 
gas accretion history and regulated by the Schmidt law. The Milky Way makes no 
exception to this behaviour. 
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1 INTRODUCTION 

The evolution of the relation between stellar mass and star 
formation rate is constrained out to high redshift (e.g., Bell 
et al., 2005, Elbaz et al., 2007, Noeske et al., 2007, Daddi et 
al., 2007, Pannella et al., 2009, Oliver et al., 2010): at any 
Z < 2, galaxies form a tight star-forming sequence, with the 
star formation rate (SFR) nearly proportional to the stellar 
mass. In particular, observations show that the average spe- 
cific star formation rate (sSFR, i.e. the SFR per unit stellar 
mass) increases by a factor of 20 from z~0 to z=2. What 
happens at z>2 is the subject of a lively debate. Whilst (ear- 
lier) studies report an almost constant sSFR (~ 2 — 3/Gyr) 
with redshift (Gonzalez et al., 2010, Reddy et al. 2012, and 
references therein), other (more recent) works support an 
increase in the sSFR with respect to the z ~ 2 value (see, 
e.g., Stark et al., 2012 and Section 3 of this article). 

The interpretation of the sSFR-redshift evolution is 
complicated by the fact that that different populations of 
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galaxies may contribute at a given redshift (e.g. Stark et al., 
2009), with Renzini (2009) arguing that the high redshift 
objects with the largest sSFR might turn into spheroids, 
whereas those with lower than average sSFR will evolve as 
discs. Indeed, a large fraction of z ~ 2 galaxies show disc-like 
structures (e.g. Cresci et al. 2009, Law et al., 2009, Gnerucci 
et al., 2011), similarly to the low-redshift counterparts in the 
mass-SFR relation. Moreover, despite these z ~ 2 galaxies 
form stars at a rate of up to a few 100 M^/yr (Shapley et al. 
2004, 2005 Daddi et al., 2007), there is increasing evidence 
that such high SFRs are not caused by merger-driven star- 
bursts as it happens locally (e.g., Rodighiero et al., 2011, 
Kaviraj et al., 2012). A large reservoir of gas, with a supply 
rate exceeding the gas consumption rate, must be available 
to these galaxies in order to maintain their high SFR (e.g. 
Erb et al., 2006, Daddi et al. 2010) for a relatively long pe- 
riod. It has been argued that, in order to do so, galaxies 
must behave as self-regulating systems that evolve in nearly 
steady-state: the gas accretion (driven by the cosmological 
accretion rate of the host halo) balances star formation and 
outflows. This is shown in several simple (analytic) models 
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of galaxy formation (e.g. Bouche et al., 2010, Papovich et 
al., 2011, Dave et al, 2012, Reddy et al., 2012, Lilly, et al., 
2013, Pipino, et al., in prep.). The above-mentioned observa- 
tions as well as the analysis performed with simple models, 
suggest that the galaxies driving the z < 2 evolution of the 
sSFR are disc-dominated. How would the Milky- Way fit in 
this scenario? What is the contribution of earlier morpholo- 
gies to the z < 2 sSFR evolution? 

If true, the constancy of the sSFR at higher redshift, 
instead, seems to require a conspiracy of many physical pro- 
cesses, including accretion, outflows, mergers and the varia- 
tion of the star formation efficiency with time and/or halo 
mass (e.g. Weinnmann et al., 2011, Khochfar & Silk, 2011, 
but see Krumholz & Dekel, 2012), to be explained in the 
Cold Dark Matter cosmology. Taken at a face value, these 
high values of the sSFR reported at z > 2, suggest a vig- 
orous star formation history, which should not differ much 
from the one inferred for massive spheroids on the basis of 
their chemical abundance pattern (e.g. Pipino & Matteucci, 
2004); however, it has yet to be proven that this is the case. 

In order to answer to these questions, we tackle the 
problem of understanding the evolution of the sSFR from the 
perspective of chemical evolution. In particular, we make use 
of models fully calibrated onz~0 observables: abundance 
ratios, ages, SN rates in the Milky Way as well as in local 
elliptical galaxies. These models well explain the observed 
evolution of the mass-metallicity-SFR relation up to z ~ 3 
(Calura et al., 2009) if one takes into account that popu- 
lations differing in morphological type contribute at differ- 
ent epochs, according to the empirical fact that the average 
age of stars in elliptical galaxies is larger than that in spi- 
rals. These models were never tuned to match the observed 
sSFR, therefore we can regard the predicted sSFR as a 
genuine result achieved without any alteration of the model 
parameters. We aim at comparing our results to available ob- 
servations, to test if the empirical trend can be reproduced 
by the star-forming progenitors of present-day galaxies of 
different morphological types. In Section 2 the models are 
briefly summarized. In Section 3 we present our results and 
discuss them. Conclusions are drawn in Section 4. 



2 THE MODEL 

The models adopted in this paper are taken from Pipino & 
Matteucci (2004) for the ellipticals, while for a Milky Way 
like galaxy we make use of the two-infall model of Chiap- 
pini et al. (1997, 2001). In brief, in all these models the gas 
accretion is regulated by an assumed infall law. Such a law 
has an exponential form for both ellipticals and spirals: 



M ac 



-t/T 



(1) 



where r is the time scale for the infall (details can be found 
in the original papers). 

Ellipticals can also eject gas (and stop star forma- 
tion) through a supernova-driven galactic wind, if the self- 
consistently derived energy budget exceeds the binding en- 
ergy of the gas. 

The star formation law is the Schmidt (1959), where 
the SFR is proportional to some power of the gas density: 



SFR 



(2) 



where v, expressed in unit of time -1 , is defined as the effi- 
ciency of star formation, namely the SFR per unit mass of 
gas, and it should not be confused with the sSFR. We adopt 
<y — o-gas (surface gas density) and the exponent k = 1.5 for 
spirals; whereas for ellipticals a is the volume gas density, 
which implies k — 1 (Kennicutt, 1998). 

The input values for relevant parameters (star forma- 
tion efficiency, infall timescale, radius) are set to those in 
Calura et al. (2009, Tables 1 and 2), where the reader can 
find also the descriptions of the star formation histories of 
these systems. In brief the Milky Way-like (massive) spi- 
ral has a central surface density Eo = 500 (2OOO)M0pc -2 , 
a scale length of 3.5 (5) kpc and a star formation efficiency 
v — 1 (2)/Gyr. The infall timescale is r ~ 8Gyr at the 
solar radius in the Milky Way spiral and increasing out- 
ward. The low (high) mass elliptical model has an effec- 
tive radius Reff ~1 (3) kpc and a star formation efficiency 
v =3 (12)/Gyr. The infall timescales are shorter (0.5-0.4 
Gyr) than in the spiral model. We adopt a Scalo (1986) 
IMF for the spirals, and a Salpeter (1955) IMF for ellipti- 
cals. However, as long as we are interested in the sSFR, 
the IMF effects basically cancel out. The present-day stellar 
masses are ~ 10 1 Mq and ~ 10 1x Mq for the low and high 
mass elliptical, respectively. The Milky- Way has a present- 
day stellar mass of ~ 10 10,2 M©, whereas the massive spiral 
has a mass of ~ 10 11 Mq. 

Chemical evolution simulations of elliptical galaxies 
have been calibrated to reproduce the observed mass- 
metallicity and mass-[a/Fe] relation in the stars of ellip- 
tical galaxies (e.g., Worthey et al 1992; Trager et al 2000; 
Thomas et al 2010, Nelan et al. 2005, Bernardi et al., 2006, 
Graves et al. 2007), under the assumption that both the 
star formation timescale and the infall timescale decrease as 
a function of galaxy mass (downsizing in star formation and 
mass assembly). These models also reproduce the observed 
chemical pattern in the gas of elliptical galaxies at both 
low (Pipino & Matteucci, 2011) and high (e.g. Matteucci 
& Pipino, 2001, Pipino et al., 2011) redshift. In particu- 
lar, Pipino et al. (2011) showed that while relatively small 
mass ellipticals well reproduce the chemical properties of 
well studied z ~ 3 Lyman Break galaxies, the most massive 
galaxies feature an abundance pattern and the dust content 
similar to those observed in the highest (z > 6) QSO hosts. 
In this paper, we model elliptical galaxies as single- zone en- 
tities of radius 10 Reff- 

As far as spiral galaxies are concerned, we adopt the 
above described Milky Way model and a model suitable for 
a more massive spiral, such as M 101, with the assumption 
that larger discs evolve faster than smaller ones (e.g. Boissier 
et al., 2001), in analogy with ellipticals and in agreement 
with observations, while still being assembled in an inside- 
out fashion (Yin et al., 2009, Marcon-Uchida & al. 2010). 
In this paper, the model predictions will always refer to the 
quantities integrated over the inner 10 kpc. 

Our models reproduce the well-established downsizing 
character of disc galaxies and spheroids, i.e. the shorter du- 
ration of the star formation in larger galaxies (Calura et al., 
09; Pipino et al., 2011). Whilst single star formation histo- 
ries predicted by the models have been shown in the original 
papers, this is the first time that we show the predicted sSFR 
for our models. 
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Figure 1. The evolution of the Milky Way (filled red circles) and 
the massive spiral model (stars and triangles) in the SFR-mass 
plane, as if they were "observed" at z=2,l and 0, respectively. 
For the high mass spiral model we show a case in which the infall 
timescale is as in Calura et al. (2009, triangles) and one in which 
it is 3 times longer (stars). The scatter in the observed relation at 
a given epoch can be explained by a small scatter in the formation 
redshift (zp) of these galaxies. In particular, for each of the two 
model for the massive spiral, we show the cases with zp = 3 
(empty symbols) and 4 (filled symbols). The relations observed 
at redshift 2 (Daddi et al., 2007, blue dashed line) and (from 
Elbaz et al., 2007, yellow dashed line) are shown for reference. 
Dotted lines give the ~0.3 dex scatter of the observed relations 
around the mean. 



3 RESULTS & DISCUSSION 

3.1 Spirals and the SFR-mass relation at z < 2 

The predicted evolution of the Milky Way model in the SFR- 
mass plane is shown in Fig. Q]by means of filled red circles. 
We also show the massive spiral model. In particular, we 
show a case in which the infall timescale is as in Calura et 
al. (2009, triangles) and one in which it is 3 times longer 
(stars). We compare the model predictions to the the rela- 
tions observed at redshift 2 (Daddi et al., 2007, blue dashed 
line) and (from Elbaz et al., 2007, yellow dashed line). 
Dotted lines give the 0.3 dex dispersion of the observed re- 
lations around the mean. The star formation history of the 
Milky Way model matches the observed evolution of star 
forming galaxies between z=2 and z=0, being almost a fac- 
tor of 2 below the average at a given mass, but still within 
the empirical dispersion. The path of a more massive spiral 
galaxy is located in the upper part of the SFR-mass plane, 
being however very similar to that of the Milky Way. For 
the generic massive spiral model we did not have any spe- 
cific constraint on the age. This allow us to show that the 
scatter in the observed relation can be explained by a small 
scatter in the formation redshift (zp) of these galaxies. In 
particular, for each of the two models for the massive spi- 
ral, we show the cases with zp = 3 (empty symbols) and 

4 (filled symbols). Our models do not incorporate mergers 
and smooth gas accretion histories are assumed, therefore we 
cannot rule out that part of the observed scatter originates 
from episodic small bursts. 

In Fig. [2] we show the temporal evolution of the sSFR. 
The light gray area brackets the observed evolution (±la) 
of the sSFR with cosmic time as compiled by Gonzalez et al. 
(2010, see also Weinmann et al. 2011) for galaxies of mass 
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Figure 2. Time evolution of the Log sSFR. The light gray area 
brackets the observed evolution (±l<r) of the sSFR with cosmic 
time as compiled by Gonzalez et al. (2010, see also Weinmann et 
al. 2011) for galaxies of mass 0.2 - 10.0 X 10 10 M Q . The larger 
dark gray box qualitatively shows the direction of the change in 
the estimated sSFR when dust extinction, correction for emission 
contamination and a large range of star formation histories are 
allowed (see references in the text and Fig. UJ. Spirals (solid: 
Milky Way; dashed: massive spiral - a model with a different infall 
timescale is drawn in green, see text for details), and ellipticals ( 
red & dashed: low mass model - red &; solid: massive model) are 
shown. 



0.2- 10.0 x lO lo Af . Let us start by discussing the z < 2-3 
evolution. At high redshift the situation is not clear-cut and 
we defer the discussion to the next session. 

Since the SFR-mass relation at z < 2 has a slope close 
to unity (e.g. Noeske et al., 2007, Elbaz et al., 2007, Daddi 
et al., 2007), the trend depicted in the figure is a reasonable 
proxy for the sSFR evolution of star forming galaxies at any 
given mass. 

We show our model spiral galaxies as black lines (solid: 
Milky Way; dashed: massive spiral). The formation redshift 
for the Milky Way model is chosen in order for the halo 
(and its associated globular cluster population) to have ages 
around ~ 12 Gyrs (De Angeli et al., 2005, Marin-Franch 
et al., 2009, Thomas et al., 2011). It is however clear that 
this model can be used as a proxy for a "normal" spiral, 
by relaxing the specific constraints on the Milky Way and 
allowing for a scatter in the values of zp as well as the infall 
timescales. For clarity purposes, we show the effect of these 
changes in the massive spiral model, where the time at which 
the star formation begins is arbitrarily chosen to be zp ~3. 
Two different infall timescales are adopted in the case of the 
massive spiral: as in Calura et al. (2009, black), and 3 times 
longer (green). 

In these models the star formation stops at the end of 
the halo phase due to the gas density dropping below the 
star formation threshold and then sets in again after the 
second infall has provided enough gas to begin the thin disc 
formation (Chiappini et al., 1997). As a consequence there is 
a short (< 0.5Gyr) period in which the star formation, and 
hence the sSFR, drops to zero. This can be seen occurring 
at about 2.5 Gyr in Fig. [2] 

The evolution of the mass-SFR relation predicted by 
our model spirals matches the observations within the scat- 
ter (c.f. Fig. [1} . Therefore it is not surprising to see that the 
temporal evolution of the sSFR at z < 2 is well reproduced. 
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More specifically, the Milky Way model evolves always be- 
low the average sSFR at any given time, but still within the 
observational 0.3 dex scatter. We do not believe that the off- 
set is significant since there is no guarantee that the Milky 
Way is the average galaxy. Furthermore, we avoided any sort 
of fine tuning. As the dashed black line shows, a later forma- 
tion redshift would automatically lead to higher sSFR that 
predicted by our best model for the Milky Way at z = 1 — 2. 
Changing the infall history (the green dashed line) makes 
the predicted values higher at z < 1. Therefore allowing for 
the scatter in these properties of otherwise similar galaxies 
can easily populate the gray area. 

The physical reason behind the agreement with data is 
that in our models, the decrease in the star formation rate 
over cosmic time is modulated by the infall law. The empir- 
ical infall law that we adopt for the Milky Way is, in turn, 
very similar to that inferred from Dark Matter only simula- 
tions, constrained to match in mass and reproduce the rel- 
atively quiet history of the Milky Way host halo (Colavitti 
et al., 2008). Such an infall law and the assumption of a 
Schmidt law for the star formation, induce the evolution 
in the gas fraction /i = a gas /atot- The predicted variation 
in the gas fraction with time is shown in Fig. O where we 
compare our spiral models to a compilation of observational 
data (Tacconi et al., 2010, Daddi et al., 2010, Geach et al., 
2010, Leroy et al., 2008). The observed galaxies have stellar 
masses above 1O 1O M0 and star formation rates such that 
they are typical star forming objects at their respective red- 
shift (i.e. they follow the observed relation shown in Fig. [TJ 
within the scatter). This ensures consistency with the other 
datasets used in this paper and with the stellar masses of 
the models. 

Since it can be easily shown that sSFR = u/x 1 /(1—jjl) 
(e.g. Reddy et al., 2012) any change in the gas fraction in 
the disc will drive a variation in the sSFR. This empirical 
finding will be demonstrated in more general and theoretical 
terms in a separate series of papers (see Lilly et al., 2013, and 
Pipino, et al., in prep.). In particular, we anticipate that a 
smoothly decreasing cosmological accretion rate, taken from 
simulations, and the regulating action of the Schmidt (1959) 
law are sufficient to drive the observed variation in the sSFR 
with cosmic time. 

Therefore we conclude that the evolution of the Milky 
Way in the mass-SFR (equivalently in the sSFR-time) plane 
inferred from models tuned to reproduce only its chemical 
evolutionary history, is consistent with the average trend of 
star forming galaxies. The predicted evolution of the sSFR 
should not be confused with the notion that the SFR in the 
Milky Way did not change by more of a factor of a few dur- 
ing its evolution, which is based on the comparison between 
the current SFR and the past average star formation rate 
(<SFR>fl 

We also note that this in agreement with earlier works 
(Calura & Matteucci 2003, Vincoletto et al., 2012), where 
we showed that the observed evolution of the luminosity 
density at z < 2 is primarily driven disc galaxies. 

Going back to Fig. [5] and the z < 2 evolution of the 
sSFR, we note that models for ellipticals are shown with 



1 Trivially sSFR ~ SFR/(<SFR> • T), where T is the cosmic 
epoch. Therefore the sSFR decreases even if SFR~<SFR>. 
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Figure 3. Gas fraction evolution predicted for the spirals (solid: 
Milky Way; dashed: massive spiral). Data on the observed molec- 
ular gas fraction from Tacconi et al. (2010, z ~ 2.5), Daddi et 
al. (2010, z ~ 1.5), Geach et al. (2011, z ~ 0.4) and Leroy et al. 
(2008, z = 0) are shown by symbols with error-bars. 



red lines. In particular, the low mass model is plotted as a 
dashed line. We show two of such galaxies, one with red- 
shift of formation of z ~ 1 (boundary of lowest ages al- 
lowed by spectral analysis of local galaxies, e.g. Thomas et 
al., 2010) and one at a higher (z ~ 3) redshift, where it 
may be observed as a Lyman Break galaxy (Pipino et al., 
2011). The predicted sSFR agrees with the observations at 
the same redshift; therefore we cannot exclude that observa- 
tional samples at z ~ 1 — 2 contain progenitors of low mass 
ellipticals. These galaxies cannot, however, explain the cos- 
mic run of the sSFR down to z ~ 0. 

The analysis of the integrated spectral light tell us that 
massive ellipticals should be already in place and/or about 
to be quenched by z ~ 2. Therefore their contribution to the 
observational samples made of z < 2 star fornring galaxies 
must be negligible. 



3.2 The specific star formation rate evolution at 

z >2 

3.2.1 Model predictions: the sSFR evolution in single 
elliptical galaxies 

On a qualitative basis, it is clear from Fig [5] that the high 
values of the sSFR, averaged over an observational sample, 
at z > 2 of ~2-3/Gyr can be easily obtained with our model 
ellipticals. In particular, the highest sSFR derived for sin- 
gle galaxies (darker gray box, see discussion below) can be 
traced back to the earliest phases of the formation of the 
most massive spheroids, whereas (lower than average) val- 
ues are typical of lower mass ellipticals, and possibly the 
halo/thick disk/bulge components of spirals. 

More specifically, our models predict rather flat star for- 
mation histories (c.f. Fig. 1 in Calura et al., 2009), irrespec- 
tive of mass, apart from a fast increase in the first 100 Myr. 
After the peak, they have a mild decline (Fig.l in Calura 
et al., 2009), certainly shallower than an exponential, un- 
til the galactic wind stops the star formation. We therefore 
expect our SFR to be qualitatively consistent with Reddy 
et al. (2012), who find that strongly declying star forma- 
tion histories give inconsistent determination of the SFR. 
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More quantitatively, the low mass elliptical model features 
SFR~ lOM^/yr, whereas for the massive one the model pre- 
dicts SFR~ lOOMe/j/r. While we cannot self-consistentljjj 
predict the SFR distribution of galaxies with our chemical 
evolution simulations, we can certainly say that these SFRs 
are reasonable in the light of those observed at z ~ 4 — 5, 
and their distribution (Smit et al., 2012). 

At the same time, the galactic stellar mass rapidly in- 
creases in the models. Therefore, for each single galaxy we 
always predict a decreasing sSFR with time; for instance 
the most massive galaxy model starts from sSFR exceeding 
10/Gyr and quickly evolves with time, much faster than in 
the case of spirals. 

This predicted behaviour is at variance with the face- 
value interpretation of the those observations (e.g. the com- 
pilation by Gonzalez et al., 2010), that seems to support a 
plateau in the sSFR. In the following sections we will show, 
however, that the plateau might arise from a combination 
of observational systematics and the combination of many 
populations of otherwise identical ellipticals that form over 
a suitable redshift range. 



3.2.2 On the observational results 

From the observational viewpoint the estimated sSFRs at 
high redshift, are subject to a number of uncertainties, and 
the average value, and the very presence of a sSFR plateau 
at z > 2 are still being debated. It is worth presenting a 
breif discussion, before making a detailed comparison to our 
models. 

In the first place, it is worth reminding that the ob- 
servational points that make the grey area in, e.g., Fig. 1, 
represents the sSFR averaged over the galaxies in each ob- 
servational sample. By inspecting the data, we derive that 
in each observational sample, there is part of the population 
with significantly higher than average sSFR, and part with 
sSFR smaller than 1/Gyr: if the proportion between these 
two populations does not change with redshift this would 
artificially create a flat trend in the average sSFR, as we 
will discuss below. 

Moreover, in those studies (e.g. Gonzalez et al. 2010) 
that use a restricted set of star formation histories in their 
SED-fitting procedure to derive masses and SFEp and no 
dust extinction, the small scatter around such mean value is 
possibly an artifact of this procedure (McLure et al., 2011). 
The debate is however still open, as it is not clear whether 
using a large range in star formation histories is physically 
sound, as, e.g. the assumption of exponentially declining his- 
tories seems to yield inconsistent results between the SFR 
estimated via SED fitting and those derived from UV and 
IR (Reddy et al., 2012). 

Additionally, an increase in the average sSFR at z ~ 
4 — 7 with respect to the z ~ 2 value (the "average sSFR" , 
lightest gray area, in Fig. [2]) has been reported by several 
very recent works: Bouwens et al. (2012), move the estimate 
of the sSFR upward by a factor of 2, due to the revised dust 



2 We would need to convolve our results with either a theoreti- 
cally motivated or an observational mass/luminosity function 

3 Which then will be correlated if, e.g., a constant star formation 
history template SED is assumed (McLure et al., 2011). 



correction. The correction for the contamination of nebular 
emission lines, that artificially increase the stellar masses 
derived from SED fitting, yield a similar up-ward increase 
in the sSFR estimate (Schaerer & de Barros, 2010, Stark et 
al., 2012). Finally, a further increase in the sSFR estimate 
is expected if exponentially-increasing star formation histo- 
rieqj are adopted (Maraston et al., 2010 but see Reddy et 
al., 2012, Gonzalez et al., 2012). 

In conclusion, the observed average sSFR at z > 4 may 
reach values of 5-10/Gyr, considerably larger than the val- 
ues observed at z ~ 2, and in better agreement with our 
predictions. 



3.2.3 Comparing models to observation 

In order to show how some of these effects would contribute 
to increase both the average observed sSFR and its scatter 
(see also Fig. 26 in De Barros et al., 2012), we first construct 
three dark gray boxes in Fig.[4l whose width along the y-axis 
is given by the la scatter around the mean sSFR estimated 
for z ~ 5 galaxies when dust extinction and a large range of 
star formation histories are allowed. In particular, we make 
use of the values published by Yabe et al. (2009), whose sam- 
ple feature a rather large average sSFR. Other observational 
samples, including those with dust correction (Bouwens et 
al., 2012) and emission line corrections (e.g. Schaerer & de 
Barros, 2010, Stark et al., 2012), would be bracketed by the 
Yabe et al. (2009) values and the Gonzalez et al. (2010) com- 
pilation (lightes gray box labelled as "average sSFR"). We 
repeat the exercise for three stellar mass bins in which we 
split the Yabe et al. data-set (different shades of dark gray, 
see Figure caption). A strong decrease in the sSFR with stel- 
lar mass becomes evident in the data. A caveat that applies 
to the data-model comparison is that the observations are 
based on UV-selected galaxies, therefore low mass galaxies 
will tend to have the highest sSFRs (e.g. Reddy et al., 2012). 

The width of the box along the abscissa is suitably cho- 
sen to bracket the z ~ 4 — 7 redshift range where, according 
to the latest data (Bouwens et al., 2012, Schaerer & de Bar- 
ros, 2010, Stark et al., 2012), the systematic uncertainty in 
the sSFR may amount to a factor of 2 or larger and the 
spread in the single galaxy measurements may be signifi- 
cantly larger than 0.3 dex. 

Let us now compare these data to the models. As a first 
step, mostly for visualization purposes, let us assume that 
all galaxies at high redshift are progenitors of massive ellip- 
ticals. This is clearly an extreme assumption, and in reality, 
as discussed below, each observational sample is made of 
galaxies of different morphologies and masses, distributed 
according to the mass function, and the observational se- 
lection functions. Therefore the average sSFR should be a 
weighted average of those predicted by different models. 
While we expect that (massive) ellipticals are important 
contributors at these redshift, a proper comparison needs a 
full self-consistent simulation of a cosmologically significant 
volunre. 

The star forming phase of the progenitors of massive 



4 That is star formation histories that are consistent with the 
reported constant sSFR 
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Figure 4. Time evolution of the Log sSFR at z >2. The light 
gray area (labeled "average sSFR") brackets the observed evolu- 
tion (±lrr) of the sSFR with cosmic time as compiled by Gon- 
zalez et al. (2010, see also Weinmann et al. 2011) for galaxies of 
mass 0.2 — 10.0 X 10 10 Mq. The three darker gray boxes show 
the sSFR estimated for z ~ 5 galaxies when dust extinction and 
a large range of star formation histories are allowed (Yabe et 
al., 2009). In particular, the darkest box (highest sSFR) corre- 
sponds to masses < 10 10 Mq, whereas the lightest of the three is 
for galaxies with masses above 1O 11 M0. The evolutionary track 
of the high mass model elliptical is rendered with a darker red in 
the phases of its evolution characterized by the stellar mass in the 
range 0.2 — 10. X 10 10 Mq. Such a model is replicated by simply 
assuming different formation redshifts zp , in such a way that the 
onset of the star formation is spaced by 0.2 Gyr. The vertical 
solid line shows how an hypothetical observation at z ~ 6 will 
read different values of the sSFR of galaxies differing only in z f ■ 



galaxies must be short (< 0.6 Gyr) to ensure that the chem- 
ical composition of their stars is strongly a-enhanced. As 
discussed above, for each single galaxy, the model predicts 
a quickly declining sSFR. The fastest decrease is for the 
most massive ellipticals. Therefore, a galaxy observed at 
say z = 6, a few Myr after the beginning of the formation, 
will have log sSFR/Gyr > 1. The same galaxy will have 
ZogsSFR/Gyr ~ 0.3 by z ~ 5 and it will be passively evolv- 
ing at z = 4. These galaxies exhibit for most of the time a 
high sSFR consistent with that observed. In particular, the 
high mass model elliptical is rendered with a darker red in 
the phases of its evolution characterized by the stellar mass 
in the range 0.2 — 10.0 x 1O 1O M and it well matches the 
data the Yabe et al. (2009) sample. 

How does this rapid decline compare to the seemingly 
constant sSFR with redhift? This might be explained by the 
fact that the formation of ellipticals of a given mass seems 
to be rather synchronized (e.g. Bower et al., 1992, Thomas 
et al., 2010). Therefore, in Fig. U we show several evolution- 
ary tracks created by assuming that the formation redshift 
of (otherwise identical) high mass models for elliptical galax- 
ies are arbitrarily spaced by 0.2 Gyr. A Szf ~ 1 Gyr scatter 
in their formation epoch is fully consistent with classic argu- 
ments based on the scatter of the color-magnitude relation 
of local ellipticals (Bower et al., 1992). The assumed redshift 
of formation are also fully consistent with spectral analysis 
of local massive ellipticals, which yields ages of at least 8 
Gyr (e.g. Thomas et al., 2010). The flat probability distri- 
bution in zf, instead, is just an assumption for illustration 
purposes, and probably does not hold in the real Universe, 
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log M, (M e ) 

Figure 5. sSFR evolution as a function of the stellar mass at 
the time of observations. Red thick lines are model ellipticals 
(solid: massive; dashed: low mass; dotted: additional case with 
~ 1O 12 M0 final stellar mass). For comparison, the massive spiral 
model is shown by a thick (dark) solid line. The black thin lines 
are the relations expected for constant star formation histories 
(500, 50, 5 M^/yr from top to bottom). Data from Yabe et al. 
(2009 - crosses), Gonzalez et al. (2010, asterisks), McLure et al. 
(2011, lozenges - converted to Salpeter (1955) IMF). 



and does not take into account the any possible diffence in 
the comoving number density of these objects at different 
redshifts. 

This exercise helps visualizing the fact that, otherwise 
similar galaxies, differing only in zp and observed at, e.g., 
z = 6 (vertical line in Fig. [4| , will display a range of sSFR 
values. That is, if one allows for a minimal scatter in the 
formation epoch of these objects, at each redshift younger 
- corresponding to higher sSFR - and older (lower sSFR) 
objects would coexist. If the probability of formation is uni- 
formly spread over the redshift range z — 4 — 10, and the 
progenitors of the massive spheroids dominate the general 
population of star forming galaxies and behave as in our 
models, then it is easy to see that at any given redshift, 
an observational sample will have the same mixture of old 
and young spheroids, yielding always an "average" sSFR 
of ~ 2 — 3/Gyr, namely the values around which our model 
galaxies spend most of their life. A series of independent gen- 
erations of galaxies has been also suggested by the compari- 
son between the evolution in the the UV luminosity function 
and the mass-UV luminosity relation of high redshift objects 
(Stark et al., 2009). 

The situation can be more complex than this, as pro- 
genitors of galaxies with different final stellar masses may 
contribute to the observations conducted at a given redshift. 
From a qualitative point of view, we can make one futher 
step and imagine that the entire population of star forming 
galaxies is made by two sub-populations: one represented 
by our low mass elliptical model, forming stars at about 
~ lOM /yr and with (relatively speaking) lower sSFR, and 
another represented by the high mass elliptical model, with 
higher (specific) star formation rates (c.f. Fig. [5] discussed 
below). Since the number densities of objects with star for- 
mation rates ~ 10MQ/yr is 0.8 dex higher than that of 
galaxies forming stars at ~ 100A/©/yr (Smit et al., 2012), 
and using the sSFR predicted by our models, we can roughly 
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estimate that the average sSFR over such a mixed popula- 
tion will be ~ 3/Gyr. 

Moreover, the decrease of the number density of highly 
star forming objects at z > 5 (Smit et al., 2012), tells us 
that we cannot assume a uniform formation probability for 
the massive ellipticals at z > 4. It may decrease with time, 
leaving us with the conclusion that at the highest redshifts 
the smaller ellipticals contribute most to the observed sSFR. 
Mergers - not included in our models - may then make some 
of these galaxies create a more massive one. Whilst it is now 
well established that late time merger-induced star forma- 
tion of (Pipino & Matteucci, 2006) and also an excess of 
dry-mergers of ellipticals (Pipino & Matteucci, 2008, Pipino 
et al., 2009), do not reproduce the abundance pattern of 
present-day ellipticals, less clear is the difference of mergers 
and revised monolithic models in the high redshift high star 
forming systems. 

The main massage of the simple exercise discussed 
above is that we showed that achieving a plateau (if any) in 
the sSFR-time evolution from the point of view of our mod- 
els is relatively easy, once reasonable constraints on zf and 
Szf are assumed even if the sSFR steeply declines with time 
for each single galaxy model. Hence it should be straightfor- 
ward for any other kind of numerical simulation, once it 
features the star formation histories for elliptical galaxies 
similar to the ones predicted by our model. The non trivial 
task is, instead, to achieve the star formation histories for el- 
lipticals with suitably high sSFR (and high [a/Fe] rations in 
the stars) from first principles (Pipino et al., 2009, Sakstein 
et al., 2011). 

Finally, our model for elliptical galaxies suggests that 
the galaxies with the highest final mass always feature the 
highest sSFR at any time along their evolution. This is 
shown in Fig. [5] where we plot the sSFR evolution as a 
function of the current stellar mass for our model ellipticals 
(red lines - solid: massive; dashed: low mass; dotted: addi- 
tional massive case built to match the highest sSFR values 
reported by observations). The green crosses are the data 
from Yabe et al. (2009, 2 = 5); whereas red asterisks are 
from Gonzalez et al. (2010, z = 7) and the purple points are 
from McLure et al. (2011, z ~ 6). The black thin lines are 
the relation expected for constant star formation histories 
(500, 50, 5 Mq/x/t from top to bottom), which are adopted 
by Yabe et al. and Gonzalez et al., to simultaneously infer 
masses and star formation rates. Model ellipticals show a 
trend which is not very different, despite having star forma- 
tion rates first increasing and then decreasing with time. 

For comparison, the massive spiral model (only disc 
phase) is shown by a thick solid line. As clear already from 
Fig. [2] the early disc phase of a spiral has similar sSFR (at 
a given mass) of those predicted by the low mass elliptical 
model. At a given mass, there is a large scatter in the data, 
with our high mass and low mass elliptical's tracks setting 
the boundaries to the observed values. 



the chemical properties of high redshift star forming objects. 
Our main results can be summarized as follows: 

• Models for spiral galaxies, where long infall timescales 
and low star formation efficiencies allow a gentle decrease 
of the gas fraction with time, reproduce the evolution in 
the SFR-mass relation from z=2 to z=0. In particular, they 
reproduce the observed factor of ~ 20 decrease in the sSFR 
since z = 2. 

• The Milky Way makes not exception to this behavior, 
evolving just above the average relation (but within the ob- 
servational scatter). 

• The halo and thick disc phase build-up implies high 
sSFR values, compatible with those observed at z > 2. 

• The early phases in models for massive ellipticals, which 
require strong burst to explain the high [a/Fe] ratios in their 
stars, explain galaxies at z > 2 with sSFR well above the 
observational average value. 

• Low mass ellipticals, whose star formation efficiencies 
(and hence sSFR) are lower than those of massive ellipticals, 
may contribute to the lower than average sSFR observed at 
z = 1 — 3 as well to the "transition" phase. Given the simi- 
larities in ages and chemical abundance pattern with ellip- 
ticals, we argue that bulges of spirals, not explicitly studied 
here, may also contribute. 

• Whilst in each single galaxy the sSFR decreases rapidly 
and the star formation stops in <1 Gyr, the combination of 
different generations of ellipticals of different mass in forma- 
tion might result in an apparent lack of strong evolution of 
the sSFR (average over a population) at high redshift. 

We thus conclude that differences in ages, chemical com- 
position and star formation timescales between spheroid- 
dominated and disc galaxies, may coherently explain the two 
phases of the sSFR evolution across cosmic time as due to 
different populations of galaxies. At z > 2 the contribution 
comes from spheroids: the progenitors of present-day mas- 
sive ellipticals as well as halos and bulges in spirals (which 
contribute with average and lower-than-average sSFR). We 
argue that a milder evolution of the sSFR at high redshift 
may be due to simultaneous observation of different gener- 
ations of ellipticals in formation. The z < 2 decrease is due 
to the slower evolution of the SFR and the gas fraction in 
discs. The same physical processes make our suite of models 
simultaneously reproduce both the mass-metallicity and the 
mass-SFR relations (Mannucci et al., 2010, Lara-Lopez et 
al., 2010) observed at different cosmic times, if we consider 
different classes of galaxies contributing at different epochs 
(proto-spheroids at high redshift, discs at z < 2, Calura et 
al., 20090 
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4 SUMMARY AND CONCLUDING REMARKS 

We show the sSFR evolution predicted by chemical evolution 
models fully calibrated on the elemental abundance pattern 
of present day ellipticals and spirals, which also reproduce 



5 This result has been recently independently confirmed by other 
models (Magrini et al., 2012). 
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